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Abstract 

Dark matter is examined within the "golden region" of the Minimal Supersymmetric Standard 
Model. This region satisfies experimental constraints, including a lower bound on the Higgs mass 
of 114 GeV, and minimizes fine-tuning of the Z boson mass. Here we impose additional constraints 
(particularly due to experimental bounds on b 57). Then we find the properties of the Dark 
Matter in this region. Neutralinos with a relic density that provides the amount of dark matter 
required by cosmological data are shown to consist of a predominant gaugino (rather than higgsino) 
fraction. In addition, the U(l)y gaugino mass parameter must satisfy Mi < 300 GeV. 
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I. INTRODUCTION 



The Minimal Standard Supersymmetric Model (MSSM) is the simplest supersymmetric 
extension beyond the standard model of particle physics, and stands to be tested in the 
upcoming Large Hadronic Collider experiments at CERN. The MSSM not only addresses 
fundamental questions in particle physics, but also naturally provides a compelling dark 
matter candidate, the lightest supersymmetric particle (LSP). In particular, the neutralino, 
which is a linear combination of the supersymmetric partners of the photon, the Z boson, 
and the neutral scalar Higgs particles, has the right cross section and mass to automatically 
provide the observed density of cold dark matter in the universe. According to the analysis 
in the latter has the value 

fi^/i^ = 0.1143 ±0.0034. (1) 

Here subscript x refers to neutralinos, h is the Hubble constant Hq in units of 100 km/s/Mpc, 
and Vt^ = p^/ pc is the fraction of the neutralino density p^ in units of the critical density 
Pc = SHq/ISttG) ~ 10~^^ h"^ g/cm^ (alternatively, Q^h!^ is the neutralino mass density in 
units of 18.79 yg/m^). 

Perelstein and Spethmann [2|] examined a particularly interesting region of MSSM param- 
eter space which they dubbed the "golden region." They argued that data and naturalness 
(i.e. a low degree of fine-tuning) point to a region within the Higgs and top sectors where 
the experimental bounds from non-observation of superpartners and the Higgs boson are 
satisfied and fine-tuning is close to the minimum possible value. They found that, in this 
region, (i) the two stop eigenstates have masses below 1 TeV, (ii) there is a significant mass 
splitting between the two stop mass eigenstates, typically 6m > 200GeV, and (iii) the stop 
mixing angle must be nonzero. They then suggested collider signatures of the golden region 
that may be found at the LHC 

In this paper, we further examine this golden region, with a particular focus on discovering 
the properties of the dark matter within it. We use the numerical package DarkSUSY 31 
to find the same golden region as [21], and apply additional constraints due to experimental 
bounds on 6 — > S7 as well. Then we look for the parameter regime inside the remaining 
golden region that also gives the right relic density of neutralino dark matter. 

Below we begin by reviewing the boundaries of the golden region, and then turn to the 
properties of the dark matter within it. 



II. THE MSSM GOLDEN REGION 



We work in the framework of the Minimal Snpersymmetric Standard Model (MSSM). For 
practical reasons, and to match the choices in [2] , we impose the following restrictions on the 
MSSM parameters: (1) we assume that all soft parameters are flavor-diagonal, (2) we assume 
a common soft mass for the first and second generation squarks, rriq = mQi,2 = m[/i,2 = 
m£,i,2, and for all sleptons, = 711^1,2,3 = m£;i,2,3, (3) we set all tridiagonal terms A to zero 
except for At, (4) we further assume that the third generation soft mass m^a = rriq, but 
let niQa and rriua vary independently. We remain with 11 free parameters: the Higgs mass 
parameter /i, the mass of the CP-odd Higgs boson, the ratio tan f3 of the Higgs vacuum 
expectafo,. values, the gaug,„o ,„ass parameters M,. M3. the soft parameters ,r, and 
mi, and the third-generation soft parameters mga, mu'i and At. Furthermore, following 
we replace the last three parameters {rnqs, mjjs and At) with the mass of the lightest stop 
mi, the mass difference between the stop masses 5m = m2 — mi, and the stop mixing angle 
6t. Finally, we use mt = 174.3 GeV for the top-quark mass, and define all parameters at the 
weak scale. 

We scan the 11-dimensional parameter space by generating random values of the pa- 
rameters. In most of our results, we use the parameters within the following ranges (all 
dimensionful parameters are in GeV): 

80</i< 500, 100 <mA< 2000, tan/3 = 10, (2) 

100 < Ml < 400, 100 < Ms < 2000, 100 < M3 < 2000, (3) 

100 < mq < 2000, 100 < m^- < 2000, (4) 

100 < mi < 1000, 100 <Sm< 600, 9t = n/A. (5) 

Notice that we fixed the value of 6t and tan/3 to reproduce one of the panels in Fig. 2 of 
2I. We define our "default scan" to be the case where these values are fixed at tan/? = 10 
and 9t = 7r/4. We also produced other special scans: we randomized tan/? in the range 0.5 
to 30; we implemented the GUT relation between Mi, M2, and M3; we separately extended 
Ml up to 2000 GeV. 

The points in Fig. [T] illustrate the golden region we obtain in our default scan. We have 
plotted the mass difference between the stop masses 6m = ffi2 — fhi as a function of the 
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FIG. 1: The Improved Golden Region: Mass difference between the stop masses 6m = rh2 — rhi 
as a function of the mass of the hghtest stop rhi. We have scanned SUSY parameter space 
with constraints imposed from the bound on the Higgs mass, fine tuning, and other experimental 
constraints including 6 ^ 57 to find an "improved" golden region; here we have fixed Of = j and 
tan/5 = 10 as our "default scan" in parameter space. For comparison, the area^ enclosed by the 
solid line shows the corresponding golden region previously found in Fig. 2 of [2]. Note that our 
lowest value of 6m = 150GeV is lower than that in the previous work as discussed in the text. 

mass of the lightest stop rhi. This plot can be directly compared with Fig. 2 of [2I, from 
which the solid triangular region is drawn. Using DarkSUSY, we have improved upon the 
previous work by applying more accurate calculations of the Higgs boson mass and of the 
b —>■ s'-f branching ratio, as described below. A variety of other experimental constraints are 
applied as well, using their implementation in DarkSUSY. In particular, LEP2 searches for 
direct production of charginos and stops constrain the chargino and stop masses, resulting 
in /i > 80 GeV and rhi ~ 90 GeV. The b ^ s^y decay rate j3] also constrains the golden 
region. Whereas noted that including this constraint was mostly beyond the scope of 
their paper, we have implemented this constraint throughout. Hence our results in Fig. [1] 
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illustrate the "improved" golden region in the presence of this additional constraint (and of 
a more accurate calculation of the Higgs boson masses). 

We note that defined a benchmark point in their Table 1 , with particular choices of the 
other parameters, and did include b —>■ s'-f ratio when they checked this point using SuSpect. 
SuSpect gave them an acceptable b ^ s'j ratio for their benchmark model, but DarkSUSY 
(which includes better expressions for the branching ratio with NLO corrections) gives an 
unacceptable value. 

The golden region in Fig. [1] contains 7000 points satisfying all relevant bounds. The 
region is a triangle, with the lower boundary due to the constraint on the Higgs mass, the 
leftmost boundary due to the bounds on the p parameter, and the upper boundary due to 
fine-tuning. We now discuss each of these bounds in turn. 



A. Lower boundary of golden triangle: Higgs mass 

The lower boundary of the golden region triangle is set by bounds on the Higgs mass. 
The LEP2 lower bound on the standard model Higgs mass is js] 

m(/i°) > 114GeV. (6) 

For generic MSSM parameter choices, the limit on the lightest Higgs is very close to this value 
as well and we may use this bound. At tree level, the MSSM predicts m{h^) < mz\cos2j3\, 
so that large loop corrections are required to satisfy this bound. The dominant one-loop 
corrections are from top and stop loops. The numerical package FeynHiggs [s] is incorporated 
into DarkSUSY to properly compute the Higgs mass and to apply the experimental bound. 
We note that the resulting boundary to the golden region that we find is slightly different 
from that of j^: e.g. at mi = 440GeV, our lowest value of 6m = ISOGeV is quite a bit 
lower than their lowest value of 6m ~ 280. The reason for this discrepancy is that their 
analytic approximations for the Higgs masses are less accurate than the values we find using 
the numerical package. 
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B. Upper boundary of golden triangle: Fine- Tuning Constraint 



Fine-tuning of the Z mass also constrains the Higgs sector. At tree level, the Z mass in 
the MSSM is given by 
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Since one of the two CP-even Higgs masses must satisfy ^ < for electroweak symmetry 
breaking, and since experimentally it is found that at least one of m^ ^, ^ m^, cancel- 
lation of the terms on the right hand side is required in order to get the right value of m^. 
Following Barbieri and Giudice j^], one may quantify this fine-tuning by computing 
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where ^ = m^, m^, b, fi are the relevant Lagrangian parameters. Then 



A{b) 



rrir. 



'l + ^^i±^tan^2/ 



1 + 



tan^ 2/3, 



u 

- cos 2/3 + ^ cos^ 13- ^ 
^ mi 



1 



^ cos 2/3 + —|- sin^ /3 - 



rrir. 



1 + 



cos 2/3 
1 

cos 2/3 



+ 
+ 



tan^ 2/3 



rriA + rriy ^ ^ 
^ \ ^ tan^ 2(3 
mi 



(10) 



where it is assumed that tan/3 > 1. The overall fine-tuning A is defined by adding the 
four y4's in quadrature; values of A far above one indicate fine-tuning. Following [2], we 
require A < 100, corresponding to fine tuning of 1% or better; this bound is implemented 
in our work and produces an upper bound fi ^ 440 GeV for tan/3 > 0.5. This matches the 
upper bound on fi found in 2|]. We impose the constraints on the chargino mass from LEP2 
chargino searches, which select ^^80 GeV. 

The upper boundary of the golden region triangle is determined by further restrictions on 
the fine-tuning. Quantum corrections to naturalness also constrain the size of the quantum 
corrections to the parameters in Eq. ([7]) . Following , here we consider the largest correction 



in the MSSM, namely the one-loop contribution to the parameter from top and stop 
loops: 

Here rrit is the top mass, A is the scale at which the logarithmic divergence is cut off, and 
finite (matching) corrections have been ignored. The correction to the Z mass induced by 
this effect is 

« (l - ^) . (12) 

To measure the fine-tuning between the bare (tree-level) and one-loop contributions 
introduced 
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Choosing the maximum allowed value of selects a region in the stop sector parameter 
space, (rrii, rh2, Ot), whose shape is approximately independent of the other parameters. This 
constraint is outlined by the upper edge in Fig. [1], which corresponds to A^ < 33.3 (3% fine 



tuning). Note that the particular va! 
TeV in this figure. As pointed out in 



ue of At depends on the scale A; we choose A = 100 
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the shape of the A^ contours and the obvious trend 



for fine tuning to increase with the two stop masses is independent of A. 
C. Left Boundary of Golden Triangle: p parameter 

The left boundary of the golden region triangle is set by measurements of the p parameter, 
which obtains corrections from stop and sbottom loops. We compute the p parameter using 
DarkSUSY and require 

(2 - 8) X 10"^ < p - 1 < (2 + 8) X 10-^ (14) 
which represents the 2a range from 81]. 

III. DARK MATTER 

Now that we have found the improved golden region with the 6 — 57 bound implemented, 
we can investigate the properties of the dark matter in this region. Using DarkSUSY, we 
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find the neutralino relic density for eacli set of MSSM parameter values in the improved 
golden region. 

Fig. [2] shows the relic density as a function of neutralino mass for all our default points in 
the golden region. As discussed previously, our default scan is defined by fixing tan/3 = 10 
and 9t = 7r/4 and scanning over other parameters. The lines illustrate the band that 
satisfies the cosmological requirement of Q^h'^ in Eq.([T]). Notice that there are points with 
the correct density for practically all neutralino masses in the golden region (with perhaps a 
tiny exception at the largest masses). Although most of the points in the figure have small 
Q-^h'^, we remind the reader that the density of points, in this and all the other plots is 
arbitrary. The simulation uses a random number generator in the parameter domain. The 
dots are used for a plot as long as they satisfy the criteria for the golden region. Thus, 
the dot density doesn't necessarily have a physical meaning. It simply represents how the 
random number generator creates the dots. 

In our default scan, we allowed the gaugino mass parameters Mi^2,3 to vary independently. 
If we instead impose the GUT relations 



we still find points satisfying the cosmological constraint on Q^h'^. This is shown in Fig. [3l 
obtained by using the parameter ranges in Eqs. [2ll5]but with the additional GUT conditions 
on Ml, M2, and M3 imposed. 

To understand the properties of the points with the cosmologically interesting values of 
Q^h"^, we have analyzed the dependence of Q^h'^ in our default sample on all of the 11 
independent parameters in the Lagrangian. Most of the parameters showed no interesting 
connection with VL-^h?, except for Mi, and the gaugino fraction Zg. 

Fig. m shows the relic density as a function of Zg/{1 — Zg) where Zg is the gaugino 
fraction of the neutralino. The denominator 1 — is the higgsino fraction. Points in the 
cosmologically interesting band typically have Zg/(1 — Zg) > 1. Thus, the dark matter in 
the golden region that satisfies Eq.(l) is predominantly gaugino rather than higgsino. 

Figs. 5(a) and (b) plot the relic density as a function of Mi, the mass of the U(l)y 
gaugino. Fig. 5(a) is for our default scan, while Fig. 5(b) is for an extended scan of 1,700 
points in which tanP varies in the range 0.5 to 30 and Mi is allowed to be as large as 2000 
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FIG. 2: Neutralino relic density ^^^^^ ^ function of neutralino mass rriy, for our default scan in 
parameter space (tan/5 = 10, 6t = -7r/4). Each dot represents a point in supersymmetric parameter 
space that lies within the golden region. The horizontal band shows the 2a range in the measured 
value of the cosmological density of cold dark matter [l!]. Notice that there are points falling into 
the cosmological band for < 300 GeV. 

GeV. One can see that points with the right relic density to explain WMAP have Mi < 300 
GeV. 

IV. CONCLUSION 

In conclusion we have imposed further experimental bounds {b 57) on the golden 
region found by [2i| and have searched for subsets of this region which provide today's dark 
matter density in the form of neutralinos. We found that Mi < 300GeV is required, and 
that the neutralino that is the dark matter is predominantly gaugino. In the future, we 
plan to examine direct and indirect detection rates in concert with LHC tests of the golden 
region. 
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FIG. 3: Same as Fig. [2] except that Mi, M2, and M3 are related through the GUT relations, 
Eqs. [T51I16I Notice that there are points falling into the cosmologically-interesting band. 
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FIG. 4: Neutralino relic density ^^h"^ as a function of j^y' '-'^^ default scan in parameter 
space. Each dot represents a point in supersymmetric parameter space that lies within the golden 
region. The horizontal band shows the 2a range in the measured value of the density of cold dark 
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FIG. 5: Neutralino relic density O^/i'^ as a function of Mi for (a) our default scan in parameter 
space, and (b) an extended scan with 0.5 < tan/? < 30 and 100 GeV< Mi < 2000 GeV. Each 
dot represents a point in supersymmetric parameter space that lies within the golden region. The 
horizontal band shows the 2a range in the measured value of the density of cold dark matter [l| . 
Notice that there are points in the golden reg^8n with the right cosmological neutralino density 
provided Mi < 300 GeV. 



